Statistical description of small quantum systems beyond weak-coupling limit 
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An explicit expression is derived for the statistical description of small quantum systems, which 
are relatively-weakly and directly coupled to only small parts of their environments. The derived 
expression has a canonical form, but is given by a renormalized self-Hamiltonian of the studied 
system, which appropriately takes into account the influence of the system-environment interaction. 
In the case that the system has a narrow spectrum and the environment is sufficiently large, the 
modification to the self-Hamiltonian usually has a mean-field feature, given by an environmental 
average of the interaction Hamiltonian. In other cases, the modification may be beyond the mean- 
field approximation. 
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I. INTRODUCTION 



03.65.Yz 



Statistical description of small quantum systems, 
which have non-very-weak interaction with environ- 
ments, is of interest in a variety of fields, ranging from the 
field of ultra-cold atoms, to chemical physics, condensed- 
matter physics, bio-physics, and so on. Recently, devel- 
opment of technologies has made it possible to perform 
direct experimental studies on thermalization of systems 
of the mesoscopic, even microscopic scale 0, H|, show- 
ing deviation from the standard canonical distribution 
due to non-weak interaction Meanwhile, theoretical 
studies and numerical simulations have been also carried 
out, some showing approach to ordinary thermal states 
0, [j] , while some suggesting significant deviations in cer- 
tain cases [B-fi"lj|. 

Although known for a long time that an equilibrium 
state of a quantum system can be described by a canon- 
ical distribution when the system-environment interac- 
tion is sufficiently weak [l2l. Il3lh a sound foundation has 
been established only recently |14| - |23j . A further ques- 
tion arises for small quantum systems, whose interaction 
with environments is usually non-negligible: In which 
way does the system-environment interaction influence 
the statistical description of such a small quantum sys- 
tem? Up to now, an explicit and general expression for 
the statistical description of such systems has not been 
derived, yet, from the microcanonical-ensemble descrip- 
tion of total systems. 

As an example, one may consider a spin at a given 
site, whose interaction with surrounding spins as the en- 
vironment is described by an interaction Hamiltonian 
H 1 = eJ s J A , where e is a parameter characterizing the 
strength of the interaction and J s and J A are operators 
acting on the Hilbert spaces of the system and the en- 
vironment, respectively The problem we are to address 
is: When e is not very small, i.e., beyond the regime of 
weak-coupling limit in which the canonical distribution 



has been derived [14Hl6l |. does the spin still have a sta- 
tistical description of the canonical form and, if yes, in 
which way does the interaction term contribute to the 
description? 

Interestingly, a similar situation is met in the field 
of decoherence when studying the so-called preferred 
(pointer) basis, in which reduced density matrices be- 
come approximately diagonal with time passing. The 
preferred basis has been known to be close to the en- 
ergy eigenbasis of the system in the case of sufficiently- 
weak system-environment interaction [23 - [2?j| . While, re- 
cently, considerable deviation has been observed due to 
non-weak interaction [27| . 

In this paper, it is shown that, taking advantage 
of techniques developed in the canonical-typicality ap- 
proach to the foundations of statistical mechanics [lj- 
|23j, an explicit expression can be derived for the statis- 
tical description of small systems, which have weak but 
not very weak interaction with environments. The de- 
rived expression still has a canonical form, but is given 
by a renormalized self-Hamiltonian of the studied system, 
which appropriately takes into account the influence of 
the system-environment interaction and reduces to the 
ordinary self-Hamiltonian in the weak-coupling limit. In 
some cases, the expression may hold even for non-weak 
interaction. 



II. MAIN RESULT 

In this section, we present our main result and discuss 
its physical meaning. (The proof will be given in the next 
section.) Before doing this, below we first introduce the 
settings and notations to be used in presenting the main 
result. 

We'll consider a total system denoted by T, which is 
composed of a small system S of interest and an environ- 
ment £ much larger than S, under a total Hamiltonian, 



H = H + H + H 



(1) 
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where H s and H £ are the Hamiltonians of S and £ de- 
termined in the weak-coupling limit and H 1 represents 
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the interaction. A parameter e is used to characterize 
the strength of the interaction, with H 1 oc e. Normalized 
eigenstates of H are denoted by \E V ) with energies E v , 
those of H s denoted by \E 8 ) with energies E^, and those 
of H £ denoted by \Ef) with energies Ef, namely, 

H\E V ) = E V \E V ), (2a) 
H S \E%) = E^E*}, (2b) 
H e \Ef)=Ef\Ef). (2c) 

We use Hse to denote a subspace of the total Hilbert 
space, which is spanned by those eigenstates \E V ) with 
energies E v lying in a narrow energy region [E, E + 8E]. 
Here, 5E is small but not too small, such that there are 
sufficiently many energy levels within [E, E + 5E] for sta- 
tistical treatment. For the simplicity in discussion, we 
assume that all the Hamiltonians discussed have non- 
degenerate spectra. 

We'll consider a renormalized self-Hamiltonian of the 
system S, as well as a renormalized interaction Hamilto- 
nian, denoted by H and H 1 , respectively, which satisfy 
the following relation, 

H s = H s + H ! s , H 1 = H 1 — H s . (3) 

Note that the total Hamiltonian H remains unchanged. 
Here, the term Hg, representing certain impact of H 1 in 
the system S, is an operator acting on the Hilbert space of 
S. We use \Es) to denote eigenstates of the renormalized 
self-Hamiltonian H s with eigenenergies E$, 

H S \E S ) = E s \Es). (4) 

(Note that, for brevity, we omit a subscript S in the 
notation |-Es).) 

The operator Hg is defined by the following relation 
in the renormalized basis, namely, 

( H a)sp = (^h, (5) 

where (H s )^ = (E 3 \H S \E^). Here, = (E^H^Ef) 
is an operator acting on the Hilbert space of the environ- 

(£) 

ment. We use T-L- to denote a subspace in the Hilbert 
space of the environment £ , which is spanned by those 
states \Ef) with energies Ef 6 [E - Ej, E - Ehy + SE]. 
The right hand side of Eq.©, namely, repre- 
sents the average of the expectation value (Ef \Hi~\Ef ) 
over \Ef) 6 H~ £ \ where 7 = /3 if E5 > Ej 3 and 7 = a if 
Ei < Ef M- 

We use p s to denote the reduced density operator of 
the system S, which is obtained from the microcanonical- 
ensemble description of the total system in the subspace 
Use, by tracing out the degrees of freedom of the envi- 
ronment. That is, 

iV5E \\E n )en SB J 



where N$e is the dimension of the subspace Use, 

N SE = £ 1. (7) 

The main result of this paper is that, under conditions 
specified below, p s has the following canonical form of 
expression with parameter /3r, 

p s ~ e~^ S /Tre'^ S . (8) 

The conditions arc: (i) The environment £ can be di- 
vided into two parts A and B, £ = A + B, such that the 
interaction Hamiltonian H 1 couples S to A only, where A 
is a small system and B is much larger than A. (ii) The 
environment £ has certain type of complexity, (iii) The 
system-environment interaction is relatively weak, (iv) 
The width SE is fixed, or decreases sufficiently slowly 
with increasing size of the environment. Details of the 
last three conditions will be given in the next section. 

To see the physical meaning of the renormalized self- 
Hamiltonian H s , let us consider the example of spin 
mentioned above. For H 1 = eJ s J A , one has (H 1 -)^ — 

eJ^~(J A )^, where j£~ = (E S \J S \E^) and (J A )^ is the 
average of the expectation value of J A in the subspace 
1-6?' . Then, making use of Eq. © , the renormalized self- 
Hamiltonian can be written as 

H s = H s + (l^J 8 + AH S , (9) 
where J A is the average of (J A )j over 7 and 

AHS = e^2(AJ A ? )J^\E 5 )(Ef\, (10) 

5/3 

with AJ A ~ = ( J A )^ - J*. The term eJ^J 5 can be inter- 

otp ' 

preted as interaction of the system S with a mean field 
J*. 

In the case that (J )j is almost the same for all the 
subspaces U ( ~ \ AJ A p ~ and H s ~ H s + eJ^"j s . For 
example, this is the case, if the environment £ is suffi- 
ciently large compared with the system S and a change 
in the energy of S has negligible influence in properties 
of £ . On the other hand, if variation of (J A )^ with 7 is 
not small, AJ A g can not be neglected; in this case, mod- 
ification to the self-Hamiltonian is beyond the mean-field 
approximation. (Generalization of the above discussions 
to the general case of H 1 = e^jJjfJi is straightfor- 
ward.) 

III. PROOF OF THE MAIN RESULT 

In this section, we give our proof for the above dis- 
cussed main result. First, we specify the second condition 
mentioned in the previous section below Eq.®, about 
complexity of £ . For the simplicity in discussion, at this 
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stage we consider a strong version of the second condi- 
tion, that is, we assume that the environment is a quan- 
tum chaotic system that can be described by the random 
matrix theory 29]. (Later, at the end of Sec lIII Al we'll 
discuss to what extent this restriction to £ can be loosed 
for the purpose here and propose a weak version of the 
second condition.) Specifically, in the expansion 

where \mA] indicate an orthonormal basis in the Hilbert 
space of the system A and \qs) for an orthonormal basis 
in that of B, and the real and imaginary parts of the 
components C l mq can be regarded as Gaussian random 
numbers with mean zero and the same variance. A re- 
lated quantity that will be used is 

<% m > = (Ef\m A )(m' A \Ef). (12) 

Next, we introduce other notations to be used in the 
proof. We use Ns, N A , Nb, and Ns to denote the di- 
mensions of the Hilbert spaces of the systems S, A, B, 
and £, respectively. We use pf E to denote the reduced 
density matrix given by a typical vector (vt^) G Use, 
namely, 

p$ E = Tr £ \*Y E )(*f E \. (13) 

Similar to H~ defined in the previous section, we use 

Ha** to denote the subspace spanned by those states \Ef) 
with energies Ef e[E-E^,E-E^ + SE] in the Hilbert 
space of the environment £ . Dimension of the subspace 

(£) (£) 

H a is denoted by N a . An important subspace to be 
considered below is Hd, defined by 

^ = ®K>^1. (14) 

a 

We use pf to denote the reduced density operator given 
by a typical vector \^> d y ) £ %di namely, 

f4=Tr e \**)(**\. (15) 

We'll consider elements of the total Hamiltonian H . In 
the basis \E^Ef), the elements are denoted by 

H ai , Pj = {E s a Ef\H\E$E^), (16) 
Similarly, for the interaction Hamiltonian, 

HiiM = (E^EflH'lElEf). (17) 
Like Hi- introduced previously, we define 

Hip = {E^\H z \El), Hl.~=(E s \H I \E£), (18) 

which are in fact operators acting on the Hilbert space of 
the small system A. Correspondingly, we have elements 



(Kph = (Ef\Hi p \Ef) and (H^ = (Ef\H^\Ef), 
and it is easy to see that 

Haifij = (H„p)ij. (19) 

The following two quantities will also be used: One is 
the following mean value, 

h% ee Tt a (H^)/N a , (20) 

the other, indicated by h, is the maximum of the mean 
value (J2mm' \ ( m A\H I a [ 3 \m A )\/ N\) for the indices a and 
P. 



A. The case of quite weak interaction 

Let us first recall two properties proved in the 
canonical-typicality approach for large environments: 

pi. p s ^ P f E M- 

P2. In the basis \E^), p% has negligibly small off-diagonal 
elements and has diagonal elements proportional to 
[l4j |. The subspaces are assumed to have dimensions 
much larger than that of the Hilbert space of the system 
S. As usually done, the environment is assumed to have a 
density of states with exponential dependence on energy 
[i~3j . then, has the usual canonical form. 

When the interaction is sufficiently weak such that the 
eigenstates lE^) are approximately equal to \E^Ef), one 
has Use — 7~td- Then, the above mentioned two prop- 
erties PI and P2 imply that p s has the usual canonical 
form. 

The above approach meets a difficulty when dealing 
with interaction with strength of practical interest, due 
to the exponential growth of the dimension of the Hilbert 
space and the at-most polynomial growth of the energy 
with increasing size of the environment, which implies ex- 
ponential denseness of the spectrum of H s + H £ and de- 
viation of T~Ld from Use- Recently, it was shown that the 
difficulty can be overcome by a perturbation approach 
and the usual canonical distribution is still obtained [l6| . 
While, it is unclear how to use this method to solve the 
previously-mentioned problem of interest in this paper 
for non-extremely weak interaction. 

Therefore, we employ an alternative method, based on 
analysis in elements of the Hamiltonians. Below, we show 
that, when the interaction is sufficiently weak (but not so 
weak such that the interaction strength is not of practical 
interest), p$ E ~ p% , then, PI and P2 imply that p s has 
approximately the usual canonical form. 

Of particular interest in our analysis are off-diagonal 
elements of the total Hamiltonian H in the basis \E^Ef), 
which are given by = (H 1 ^)^ [Eq.©]. After 

some analysis (see appendix [SJ, we find that (H^g)ij 
have quite different properties depending on whether i = 
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j or not, specifically, 

(H^) u ~ h d a % (21a) 

(^)« ~ N 3 f/N^ 2 fori^j, (21b) 

1(^)^1 < hN*f/Nl' 2 iovi^j, (21c) 

These relations imply that, in the general case with non- 
zero h^p, one has 

|(^)«| >|(^k| of i^j, (22) 

due to the largeness of the environment. Below, we dis- 
cuss this general case of non-zero h^S [30| . 

For sufficiently small e, the estimate (|21al) shows that 
IC^afl)**! can De mucn smaller than the corresponding 
unperturbed-level spacings, — E^ (the system S hav- 
ing a non-degenerate spectrum). In this case, the influ- 
ence of the elements o i in the wave function of \E V ) 
in the basis \E^Ef) can be neglected. On the other hand, 
nearest-level spacings are on average approximately pro- 
portional to 1/Ng, while (H^p)ij ~ ^/N^ 2 for i ^ j 
[see (|21bl) ]. hence, as long as the environment is suffi- 
ciently large, couplings between nearest unperturbed lev- 
els are usually much larger than their spacings. There- 
fore, for whatever small but fixed e, each \E V ) has signifi- 
cant expansions in many basis states \E^Ef), as a result, 
Hse — Ha can not hold. 

To go further, we need to get an estimate to the width 
of the wave function of \E n ) in the basis \E^Ef). As 
discussed above, in this weak-coupling case, the terms 
{H^pju can be neglected. To get the estimate, we em- 
ploy a first-order perturbation-theory treatment. The 
basic idea is that perturbation theory is applicable to 
components in those basis states whose unperturbed en- 
ergies are sufficiently far from the eigenenergy E v |3lll3^ . 
hence, perturbatively computing population of \E V ) in all 
far- lying unperturbed basis \E^Ef), one may find an up- 
per border for the width of the main body of the wave 
function. The result is that the main body of the wave 
function lies within a narrow energy window around E n , 
with width de oc 1/AE (see appendix|B|. Here, AE is the 
scope of the eigenenergies of the total system, such that 
the averaged density of states is around NsNg/AE. In 
the large-environment limit, AE goes to infinity, hence, 
the width de shrinks to zero. 

Now, we specify the exact content of the condition (iv) 
discussed below Eq.©, that is, we assume that the width 
8E is fixed, or decreases sufficiently slowly with increas- 
ing size of the environment, such that Se/SE goes to zero 
in the large-environment limit. As shown in Appendix [Cl 
under this condition, pf E has properties similar to those 
of discussed above in P2, hence, has the usual canon- 
ical form. Then, according to PI, p s also has the usual 
canonical form. 

Finally, we observe that the above-used strong version 
of the second condition, i.e., the assumption of £ being 
a quantum chaotic system, can in fact be loosed for the 



purpose of getting the above-discussed estimates to the 
elements (H^s)ij- Indeed, as shown in App[D] similar 
estimates to (H^g)ij can be gotten under the following 
assumption about the environment, namely, 

• £"mm °f a given \Ef) does not change much with m 
and | G^ m , | with j ^ i does not change much with 
m, m! , i and j. 

Obviously, this assumption is valid for the above- 
considered chaotic environment, hence, we call this as- 
sumption a weak version of the second condition. Then, 
following the same arguments as given above, one also 
finds a canonical form of p . 

Furthermore, it may happen that the above-discussed 
weak version of the second condition is satisfied only 
within a part of the spectrum of an environment of in- 
terest. For example, if the environment has a sufficiently 
chaotic motion in the middle energy region, while has a 
regular motion with integrable feature in the low energy 
region, then, the above-used approach is applicable only 
within the middle energy region. It is straightforward to 
extend discussions given in ApplDlto this case, where the 
relevant change will be that Ng in the relation (1D6[) . so 
in relations (I21b[) and (|21cp . is replaced by the number 
of the levels in the part of the spectrum of interest. 



B. The case of relatively weak interaction 

We use e ew to indicate a border of e, beyond which 
some elements (-ffig)ii become non- negligible compared 
with the corresponding level spacings. The case of e < 
e ew having been discussed in the previous section, in this 
section we discuss the case of e > e cw . In this case, argu- 
ments given in AppjB]and[Cl to show the canonical form 
of p$ E in the basis \E^), do not apply. In fact, here, a 
total state {E^) may have comparable populations in re- 
lated basis states \E^Ef) and \E§Ef), as a consequence, 
the width of the main body of its wave function does not 
shrink to zero in the large- environment limit. 

A main observation is that the renormalization of the 
self and interaction Hamiltonians given in Eq. (J3|) can sig- 
nificantly reduce the elements (HL~)u in the renormal- 

ized basis \Es). When the reduction is sufficiently large, 
arguments given in ApplBl and [Cl will become applicable 
in the renormalized basis. To show this, below we dis- 
cuss two situations separately, related to validity of the 
so-called 'eigenstate thermalization hypothesis' (ETH) 
[291 133 . The ETH, confirmed by direct numerical simula- 
tionsjl, H| , effectively states that a few-body observable 
(e.g., Hi- here) has almost the same expectation values 

for energy eigenstates within an appropriate energy win- 
dow of a large, interacting many-body system. We use 
^eth to indicate a subspace spanned by eigenstates \Ef) 
lying in an energy window of the environment, for which 
the ETH is applicable [H. 
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Let us first discuss the case that there exists a sub- 
space %eth > which is sufficiently large to include all the 
subspaces W~' defined in SecHU In this case, (HL~)ty — 
(HL~)eth, where )eth indicates the average of the 

* cap' x cap' 

expectation value (Ef\H^~\Ef) over \Ef) 6 %eth- Us- 
ing this property and the definition of Hi in Eq.©, we 
get the following explicit expression of H s , 



H S — /] (Hi^)-ET-g\Es) {Efj\ - (If 7 ) ETH, 
5/9 



(23) 



where (-ff^ETH indicates the average of (Ef\H I \Ef) over Pse 



\Ef) 6 %eth> which is an operator acting on the Hilbert 
space of the system S. Substituting the expression of H 1 
in Eq.© into the definition {H 1 -)^ = (E%Ef \ H 1 1 E^Ef ) 
and making use of Eq. (l23t and the ETH, we find almost 
zero (Hl~)u, namely, 



(H^) u = {Ef\Hi~\Ef) - (EslH'slEj;) ~ (24) 

for all \e[ £) ) e W^th- 

Due to the smallness of the elements (Hl~)a in 

Eq.(|24p. arguments used in AppjB] and [C] are applica- 
ble for the renormalized Hamiltonians and in the renor- 
malized basis. Hence, the system S has approximately 
the canonical statistical description in Eq.© in terms 
of the renormalized self-Hamiltonian. Note that, since 
Eq.(f24f holds even for not small e, in this case of suf- 
ficiently large %etH' Eq.JSJ may hold even when the 
system-environment interaction is not weak. 

(£) 

Next, we consider the case that no subspace % E th 

(£) 

is sufficiently large to include all the subspaces H~ . 
Since we have assumed validity of the microcanonical- 
ensemble description for the total system within a nar- 
row energy window 5E, it is reasonable to assume that 

(£) 

the ETH is applicable within each subspace "H- with 
the same width of energy window (the system S being 
small). Then, for E% > E^, making use of Eq.© and 

the definition of H 1 , we get 







fOT lEf) G U f (25) 

K^-K^h fOT \ E f)^f ' 



Since the dimension of is larger than that of H~ 

for E~ > E~, Eg. ([23)1 shows that most of the elements 

(HL~)ii are negligibly small. Meanwhile, since (H^)u ~ 

{H* p )y for \Ef) e r6y ] if the ETH is applicable within 

T~6y\ in general, the rest of the elements {HL~)u (the 

second line on the right hand side of Eq. ([2~5|) ') are also 
reduced compared with (H^,a)u. It is easy to see that 



,(£) 



similar results also hold for E~ < E~. To summarize, 
compared with (H^o)u in the unrenormalized case, the 
elements (H~z)u are considerably reduced. 

As discussed in the previous section, for e < e ow , 
(H~as)u are small enough to guarantee validity of the ar- 
guments used in ApplBlandlCl which predicts a canonical 
form of pf E . The above-discussed reduction of (Hz~)u 
implies that, at least for e not too much above e ow , 
(HL~)h should be small enough to guarantee validity of 

the arguments used in App IB1 and [Cl for the renormalized 
Hamiltonians and in the renormalized basis. In this case, 
g E has a canonical form in terms of the renormalized 
self-Hamiltonian, hence, Eq.© holds. This completes 
our proof. 



IV. CONCLUSION AND DISCUSSIONS 

In this paper, we have studied statistical description of 
a small quantum system S in relatively weak interaction 
with a huge environment £ . It is shown that the system 
S has a canonical statistical description given by a renor- 
malized self-Hamiltonian, which appropriately takes into 
account the influence of the system-environment interac- 
tion, if the following four conditions are satisfied: 
(i) The environment £ can be divided into two parts A 
and B, such that the interaction Hamiltonian H 1 cou- 
ples S to A only, where A is a small system and B is 
much larger than A. (ii) Eigenstates of the environment £ 
within a sufficiently large region of its spectrum have the 
property that G^ m of a given \Ef) does not change much 
with m and \G^ m , \ of j ^ i does not change much with 
to, to', i and j. (iii) The system-environment interaction 

is relatively weak, in the case that the ETH is applica- 

(£) 

ble within each subspace ~H S separately; while, there is 
no restriction to the interaction strength, if the ETH is 
applicable within the direct sum of the subspaces ~H~ . 
(iv) The width 6E is fixed, or decreases sufficiently slowly 
such that 5e/SE goes to zero in the large- environment 
limit. 

Some remarks concerning the above conditions may 
be helpful: (a) Concerning the above first condition, the 
division of £ into A and B is not a too strict restric- 
tion in interaction, since it is usually satisfied in physical 
models of practical interest. While, the first condition 
may be unsatisfied, if the system-environment interac- 
tion has a long-range feature, (b) The above second con- 
dition is satisfied by an environment which is a quan- 
tum chaotic system that can be described by the random 
matrix theory. However, the exact relation between the 
above-discussed property of G^ m , stated in the second 
condition and integrability/non-integrability of the envi- 
ronment needs further investigation, because many-body 
effects may also play important role here, (c) In the case 
that the ETH is completely inapplicable for the environ- 
ment, discussions given in Sec lIIIBl do not apply, hence, 
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it is not clear whether Eq.Q is valid for relatively weak 
interaction in this case. 

It may be also helpful to give some further discussions 
for the example of spins mentioned in the section of in- 
troduction. (1) Models like Ising model with nearest- 
neighbor coupling satisfy the above first condition about 
separation of £ into A and B\ here, A is composed of 
spins that are the nearest neighbors of the spin S. Dis- 
cussions given in this paper are not valid for the case 
that the spin S is directly coupled to all the spins in its 
environment. (2) In a model like Ising model, when the 
coupling among the spins within the environment are not 
weak, it is usually reasonable to expect that the above 
second condition may be satisfied. 

Before concluding this paper, we give a brief discus- 
sion for possible experimental test in effects of some non- 
mean-field part of the renormalized Hamiltonian, e.g., 
the term AH S in Eq.© in the case of H 1 = eJ s J A . 
The term AH S may cause notable effects, when the av- 
eraged expectation values of J A in different energy re- 
gions [E — E a , E ~ E a + SE] have notable difference. For 
example, this may happen for an environment lying in a 
state in the vicinity of a (quantum) phase transition. In 
particular, if the environment has localized states and ex- 
tended states separated by certain energy level, then, J A 
may have different expectation values below and above 
the energy level, which may have experimentally testable 
effects. Another possibility lies in that, for a relatively 
small environment (still much larger than the system S), 
it may be relatively easy for the averaged expectation 
value of J A to have notable variation with the energy 

Ea- 
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where G^ m , is defined in Eq. tTlIZ]) . Using the expansion 
in Eq. iriTj) . we get 

Gmm> = (Cmq) ^m'q- (A2) 

Normalization of \Ef) requires that 

£|Ct/ = l. (A3) 

mq 

As assumed, the real and imaginary parts of C l mq can 
be regarded as random numbers with mean zero and the 
same variance, hence, Eq. (|A3p implies that most of the 

1/2 

values of \C % mq \ are around N £ . Then, noticing the 
largeness of Nb and that Ng = NaNb, making use of 
Eq. (|A2[) and the randomness of , we get 

\Glm> I - (NsNa)- 1 / 2 for i^joTm^m'. (A4) 

For i = j and m = m! , we have 

G* m * 1/N A . (A5) 

Substituting the above-obtained estimates to G into 
Eq. (|Al[) and using the approximation (1 — N^ 2 N^ 1 ^ 2 ) ~ 
1, we get 

(Hi p ) u ~ h% + Ah, (A6) 

where h^p is defined in Eq.([2"0"|) and Ah satisfies \Ah\ < 
KpN A /2 Ng 1/2 . Here, 

(A7) 

A mm' 

i.e., the average of the absolute values of the elements of 
i?£g. Similarly, we find 

Kff^kl < h aP N 3 A /2 N- 1/2 < hN 3 A /2 Nz 1/2 for i + j, 

where h, defined below Eq. (|20p . is just the maximum of 
h a p. 



Appendix A: Estimates to (H a a)ij (I) 

In this appendix, we give estimates to the elements 
(H^p)ij = (Ef\H^JEj), under the assumption that the 
environment can be regarded as a quantum chaotic sys- 
tem that can be described by the random matrix theory. 
Since H 1 couples the system S to A only, H^p is in fact 
an operator acting on the Hilbert space of the system A 
and its elements can be written as 

(H a ph = J2( mA \ H U m 'A)(Z m >, (Al) 

mm' 



Appendix B: Width of \E V ) in basis \E%) for 
sufficiently small e 

For brevity, we use \E%) to denote the basis states 
\E^Ef) (E% = E% + Ef) with increasing energy order, 
and use C v k to denote the components of an eigenstate 
|^) in the basis |J5jg), namely, C vk = {E%\E V ). The 
main body of the wave function of \E n ) in the basis \El) 
lies within a region denoted by [k\, k?\, outside which the 
total population of \E^) is equal to some small positive 
number e p , namely, Pkxk 2 — e pj where 

P klk2 = E \C vk \ 2 . (Bl) 

keQ 
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Here and after, Q denotes the set of the label k lying 
outside the region [ki,k 2 ], i.e., Q = {k : k < k% or k > 
k 2 }. 

Below in this appendix, using a first-order perturbation 
theory and results of App[S] we derive an estimate to the 
width Se = E® — , which gives an upper border for 
the width of the main body of the wave function of \E V ). 
Using \E% ) = \E^ Q Ef Q ) to denote the basis state whose 
energy is the closest to the total energy E v , the first- 
order perturbation theory gives the following prediction 
for the population Pk 1 k 2 , which we denote by P^ k2 , 



P (i) _ 



E 

k£Q 



El -El 



with \E° k ) = \ElE £ J ). (B2) 



Note that fco lies in the region [fci , k 2 ] , hence, not belongs 
to Q. 

Then, substituting the relation (|A9[) for j ^ io into 
Eq.((B2|, we get 



p(D < 

kik 2 ~ ]\f £ 



keQ 



1 



El -El 

fc k 



Pd, 



(B3) 



where the prime means that the summation is over those 
k G Q for which j ^ io and pd is the contribution from 
k S Q with j = io, namely, 



Pd 



E 

k£Q&zj=i 



£2 



(B4) 



Noticing that the label k represents the pair ((3,j) and 
the summation over "fc S Q&j = io" is usually given by 
a summation over some of the labels f3 ^ ao, we get 



be taken away. Then, replacing the summation by an in- 
tegration over E°, namely, by J dE°p(E°), where p(E°) 
is the density of states, we get 



fclfc2 ~ 



h 2 N\ 
N e 



dE°- 



P(E°) 



E° k 



(B8) 



where the prime indicates the integration domain 
(-00, U[<, oo)^ 

To have a rough estimate, we make a further approx- 
imation that p(E°) in (|B8j) can be approximated by its 
average value within an energy region with width AE, 
which includes most of the levels of the total system, 
that is, p(E°) ~ N £ N S /AE. Then, performing the inte- 
gration on the right hand side of (|B8[) . we have 

(1) h*N A N s ( 1 , 



r kik 2 



AE 



\E° k0 -El\ 



\E° h -El 
i fc fc 2 



(B9) 

Replacing P k l\ 2 by e p , the relation (|B9[) gives the follow- 



ing estimate to 5e 



K2 



5e (1) < 



4h 2 N A N s 



e p A£: 



(BIO) 



In the limit that the B part of the environment becomes 
infinitely large, AE goes to infinity, while, both N$ and 
Na remain finite. Therefore, for each fixed e p , 8e^ — > 0, 
that is, the energy width oi\E v ) in the basis \Ei) shrinks 
to zero, in the limit of large environment. 

Finally, we give a remark on the validity of using per- 
turbation theory to give the above estimates. In fact, as 
shown in Refs.[32|, l35j . for E% and E^ sufficiently far 
away from E^, the components C n k of k € Q can be ex- 
pressed in convergent perturbation expansions in terms 
of C vk of k e [h, k 2 ] . 



Pd Pa i j 



where 



Pa i 



E 

/3(#« ) 



(Hi 



E S - E S 



(B5) 



(B6) 



Using the estimate in (|A6[) and neglecting Ah, we find 



< 



[h d fq a with q a = Yl \ E * 



E p\ ■ 



(B7) 



Here, hd is the maximum of 

When e is sufficiently small, since hd is proportional to 
e and q a is independent of e, one has p a i -C e p for a given 



e„. Then, due to the relation (IB5|) . the term pd on the 



right hand side of (|B3[) can be neglected. 

To go further, since the number of k for which j = io 
is much smaller than that for which j =^ io, the prime 
over the summation on the right hand side of (|F33|) can 



Appendix C: Canonical form of pg E for sufficiently 
small Se 



In this section of appendix, we show that pf E has the 
usual canonical form, if SE decreases sufficiently slowly 
with increasing size of the environment, such that de/SE 
goes to zero in the limit of large environment [cf. the 
inequality (|B10|) ]. 

Normalized typical vectors in the subspace Use can be 
written as 



\K E ) 



ty 



E C n \E v ), 



(CI) 



where T$e is the set of the label r\ for which \E n ) £ 
Hse, namely, Tse = {ij : \E n ) e Use}, and the real 
and imaginary parts of C v are independent real Gaussian 
random variables with mean zero and variance 1/2. Let 
us expand the eigenstates \E V ) in the basis \E^Ef), 



\E v ) = j2aE s a Ef). 



(C2) 
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Substituting Ea.([C2]l into Ea.([UT|). we write |*£f ) in a 
form similar to that of a typical vector in the subspace 
7id, which has been used in Ref.[14| in deriving properties 
of pf , namely, 



(C3) 



where 



i^hE a -i^>> 



K ai - E Cjj/ai- 

ver SE 



Noticing that K. 2 SE = J2 v er s 
forward to verify that 



|C„| 2 



(C4) 
(C5) 

Nse, it is straight- 



<«£l«S> " 



Let us first discuss the case of a — (3, in which 



(C6) 
(C7) 



(«=EE E C n C n'U V J*ft (C8) 



The main contribution to the right hand side of Eq. (|C8 
comes from the diagonal terms with 77 = rj' , hence, 



E E 

i n^TsE 



\c v \ 2 \f, 



2| f ?) |2 



(C9) 



On the right hand side of Eq. (|C9l) . with a fixed, we need 
to consider only those indices i, for each of which the 
basis state \E^Ef) lies within the main body of some 
state \E V ) of -q e r^. That is, we need to consider only 
those i, for which the corresponding energies Ef lie in 
the energy region [E - E% - Se, E + SE - E% + Se]. We 

use N^g e to denote the number of levels Ef lying in this 
energy region. 

Now, we are ready to compute {pf E ) a a- Since Se <C 
SE, for most of the indices i mentioned in the previous 
paragraph, we have 



If, 



'/ 2 



1. 



which implies that 



Another result of Se <C SE is that N 



1. 



a.Se 



(CIO) 



(Cll) 



Ni £) . Then, 



Eq.dCaJ gives (fi£|fig) « M £) , hence, 
(pf B ) QQ « N^/N SE . 



(C12) 



Next, we discuss the case of a ^ (3. It proves useful to 
explicitly express some properties related to the narrow- 
ness of the wave functions of \E V ). We use n m to denote 



the average number of major components of \E V ) in 
the basis \E^Ef). Thus, on average, |/^| » (ri m )~ 1/2 . 
The major components lie within the main body of the 
corresponding wave function with width Se. In the limit 
of large environment, since Se/ SE goes to zero, n m /N$E 
also goes to zero. 

Note that each state \E^Ef) also has on average n m 
major components, when \E n ) is taken as the basis. 
Hence, the summation on the right hand side of Eq. (|C5[) 
is effectively taken over n m major terms. Then, due to 
the randomness of the coefficients C v and the relation 
I /at I ra (^m) -1 / 2 , we have \K ai \ ~ 1. Furthermore, as 
discussed above, for a given a, the number of the indices 



{£) 

Se 



hence, 

(£) 



i needing consideration is about N£ 

the number of non-negligible K a i is also about N& 
Meanwhile, the number of the random coefficients C v is 

Nse = J2 a ■ Therefore, Kpi and K ai with j3 ^ a in 
the summation in Eq. (|G7|I can be treated as independent 
random numbers. Then, from Eqs. (|C6[) and (|CJ7|) . we get 

for a//}, 

\(p!eU\<(N se )- 1/2 . (C13) 

Equations (|C12I) and (|C13[) show that pf^ has approx- 
imately the same properties as pf stated in P2 (see 
Sec lIII Al) . hence, pf E also has the usual canonical form. 



Appendix D: Estimates to (i?„a)y (II) 

In this appendix, we show that estimates, similar to 
those given in AppJX] for the elements (Hip)ij, can be 
obtained under the assumption that G^ m of a given \Ef) 
does not change much with m and |G^ m , | of j / i does 
not change much with m, m', i and j. This assumption 
is weaker than the assumption that the environment is 
a quantum chaotic system that can be described by the 
random matrix theory. 

Let us first discuss the case of i = j. A major con- 
tribution to {H^p)a from the right hand side of Eq. dAip 
comes from the diagonal terms with m — m' , for which 
— \\( m A\Ef)\\ 2 . Since G^ m does not change much 
with m and 



E — (Ef\Ef) — 1, 



(Dl) 



we have G^ m ~ 1/Na- This gives the following estimate 
to the above-mentioned major contribution to (-ff£g)ii, 
namely, 



{m A \HL\m A )G l rt 



' dia 



(D2) 



Furthermore, since \G^ m , | with m ^ m! is usually not 
larger than G^ m , Eq. (|Al[) and the above result G" m ~ 
1/Na give the following estimate for an upper border of 
\( H ip)u\> namely, 



\{HL) vl \ < N A h. 



(D3) 



To get an estimate to |G^ m , | with j ^ i, we note that 

= (Ef\m A )(m A \Ef) ~ ^(m) 

where the identity J2j \Ef)(Ef\ = 1 nas been used. As 
assumed, |G^ m , | of j ^ i does not change much with m, 
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m', i and j, hence, Eq. (|D4j) implies that, in most cases, 
\G l l m ,\^(N £ N A )-^ iovi^j. (D5) 
Using this estimate and Eq. (|Alj) . we get 

\(H^h\ <hN 3 A /2 N~ 1/2 fori^j. (D6) 
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